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ABSTRACT 

A  linearized,  input-output  model  matching  problem  for  single-input,  multi¬ 
output  nonlinear  systems  is  formulated  in  the  frequency  domain,  and  solved 
using  dynamic  output  feedback.  The  objective  is  to  obtain  a  desired  transfer 
function  for  the  family  of  closed-loop  linearized  systems. 


Introduction 

K  In  recent  years  there  has  been  rapid  development  of  the  closely-related  pseudo- 
linearization^?^ and  extended  linearization  [1,5]  approachs  for  design  of  nonlinear  control 

.  J  - - M 

laws  for  nonlinear  systems.  In  this  paper?  we  adopt  a  frequency-domain  viewpoint,  and  use 

results  in  -{§}  to  formulate  and  solve  a  linearized  model  matching  problem  for  single-input, 

multi-output  nonlinear  systems.  Die  objective  is  to  construct  a  nonlinear,  dynamic,  output 

feedback  control  law  so  that  the  resulting  closed-loop  system,  when  linearized  about  its  family 

of  constant  operating  points,  has  a  linearization  family  with  transfer  function  that  matches  — — 
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exactly  a  given  parameterized  transfer  function.  If  the  given  transfer  function  is  parameter 
independent,  then  this  problem  can  be  viewed  as  a  type  of  input-output  pseudo-linearization 

a 

problem.  -  - _ _ 

We  consider  a  single-input,  multi-output  nonlinear  system  with  input  u(t)eR  and  out¬ 
put  y  (t)  €  Rp,  defined  for  t  >  0.  The  system  is  assumed  to  have  a  constant  operating  point 
family  { [u(a),  y(o)],  aeT}  where  the  parameter  set  TcR  is  an  open  neighborhood  of  a 
fixed  cq€R.  That  is,  for  a  constant  input  u  (/)  =  u(c*i),  where  aj  6  T  is  fixed,  there  exists  an 
initial  condition  for  the  nonlinear  system  such  that  the  output  is  the  constant,  y(t)  =  y(at). 
Typically,  the  constant  operating  point  family  is  parameterized  by  constant  values  of  the  input, 
one  of  the  outputs,  or  one  of  the  underlying  states  of  the  system. 

We  will  adopt  a  local  viewpoint  throughout  this  paper  and  implicitly  allow  shrinking  the 
size  of  r  in  order  to  economically  state  results.  (When  particular  examples  are  considered,  it 
often  is  possible  to  solve  the  problem  in  a  nonlocal  fashion.)  To  avoid  counting  the  order  of 
continuous  differentiability,  we  will  assume  that  all  functions  and  their  derivatives  that  appear 
in  the  sequel  are  continuous. 

Corresponding  to  the  operating  point  family  {[u(a),  y(a)],  oeT},  let  the  px  1 
parameterized  transfer  function  of  the  linearization  family  for  the  given  nonlinear  system  be 


-  -^4 

d(s,  a) 


(1) 


where  d  (s,  a)  and  each  entry  of  N(s,  a)  are  polynomials  in  s  with  coefficients  that  are  real 
functions  of  a,  and  d (s,  a)  is  monic.  We  assume  that  for  each  aeT(l)  is  irreducible.  Also 
we  assume  d( 0,  a)  g*  0,  and  N( 0,  a)  ^0  for  each  qGT.  Thus  Ta(s)  has  no  fixed  poles  at 
s  =  0,  and  not  all  entries  of  Ta(s)  have  a  zero  at  s  =0.  Finally  we  assume  that 


£(a)  =  ra(0)^(a),  aer 


(2) 


It  should  be  noted  that  that  if  d  (0,  a)  f  0  for  each  a  €  T,  then  (2)  is  a  sufficient  condition  for 
Ta(s)  with  {[u(a),  y(a)],  aeT}  to  be  a  linearization  family  for  a  nonlinear  system.  If  the 
linearized  state  equation  for  the  nonlinear  system  has  no  eigenvalues  at  zero  for  each  a  €  T, 
then  (2)  also  is  necessary.  [5] 


-3- 


Linearized  Model  Matching  j 

i 

The  linearized  model  matching  problem  involves  computing  a  nonlinear  dynamic  output  I 

•  i 

feedback  control  law,  with  new  input  w(t)eR,  such  that  when  the  closed-loop  nonlinear  sys¬ 
tem  is  linearized  about  its  constant  operating  point  family,  the  resulting  linearization  family 
has  a  given  parameterized  transfer  function.  The  given  parameterized  transfer  function 
(model)  for  the  linearized  closed-loop  system  is  written  in  a  form  similar  to  (1), 


JT(0- 


Nm(s,  a) 
dm(s,  a) 


I 

(3)  i 


where  the  superscript  m  stands  for  model  and  dm(s,  a)  is  monic.  We  assume  that  for  oe  r  ; 

^"*(0,  a)  fO  and  Nm( 0,  a)  j*0,  so  that  7J(j)  has  no  pole  at  s  -  0  and  not  all  entries  of 
T£  (s,  a)  have  a  zero  at  s  «  0.  Also  we  assume  that  a  family  w(a)  of  constant  operating  point 
values  for  the  new  input  is  given  such  that 

i 

£(a)»7?(0)^(a),  a€T  .  (4)  i 

so  that  T£(s)  with  {[w(a),  y(a)],  a  €  T}  is  a  linearization  family. 

The  linearized  model  matching  problem  will  be  solved  as  follows.  First  we  compute,  if 
possible,  a  parameterized  scalar  linear  precompensator  T%(s)  such  that  ! 

Ta(s)ra(s)  =  7 ?(s)  (5) 

This  is  illustrated  in  Figure  1  where  the  subscript  S  indicates  deviations  from  constant- 
operating-point  values.  Then  we  implement  the  precompensator  T%(s)  using  a  parameterized 
linear  dynamic  output  feedback  configuration.  Finally,  a  nonlinear,  dynamic,  output  feedback 
law  is  computed  whose  family  of  linearizations  is  precisely  the  parameterized,  linear,  dynamic, 
output  feedback  law. 

From  (5)  it  is  clear  that  N(s,  a)  and  Nm(s,  a)  must  satisfy 

b(s,  a)N(s,  a)  =  a(s,  a)Nm (s,  a),  qGT  (6) 

for  a  (s,  a)  and  b  (s,  a)  which  are  scalar  polynomials  in  s  with  coefficients  that  are  real  func¬ 
tions  of  a.  If  (6)  is  satisfied,  the  parameterized  linear  precompensator  will  have  the  form 


7 x  b  (s,  O t)  d (s,  Or) 

aU  a(s,a)dm(s,a)  ’ 


ae  r 


(7) 


For.7£ (s)  to  be  proper  we  must  have,  for  each  aeT, 


6d  (s,  a)  -  6iN(s,  a)  <  6dm(s,  a)  -  6iNm(s,  <a) ,  i  *  1, ....  p  (8) 

where  S  denotes  polynomial  degree  and  Si  denotes  polynomial  degree  of  the  t*  row  of  a  poly¬ 
nomial  column  vector.  For  example,  if  the  plant  and  the  model  are  such  that 
Sd(s,  a)  *  6dm(s,  a),  then  (8)  implies  that  there  should  be  no  more  finite  zeros  in  the  model 
T%  (s)  than  in  the  linearized  plant  Ta(s). 

The  precompensator  7^  (s)  is  to  be  implemented  by  a  parameterized,  linear,  dynamic, 
output  feedback  law.  Although  different  implementations  are  available,  a  parameterized  ver- 

t 

sion  of  the  implementation  given  in  [2,  pp  519-522]  will  be  developed  This  is  presented  in 
some  detail  because  of  additional  hypothesis  required  due  to  parameterization  by  a.  < 

Step  1:  Assuming  that  Sa  (s,  a)  in  (6)  is  constant  for  a  €  T,  compute  least  degree  polynomials 
rtf(s,  a)  and  df(s,  a)  such  that 


rtfjs,  a) 
d/(s,  a) 


Tjjs)  =  b(s,a) 
d(s,  a)  a(s,  a)dm(s,  a)  ’ 


aeT 


(9) 


Step  2:  Letting  6dj{s,  a)  =  /  and  6d(s,  a)  =  n,  we  set 

df(s,  a)  =  v(j,  a)df{s,  a)  ,  Ufa,  a)  =  v(s,  a)n^s,  a)  (10) 

where  v  (j,  a)  =  1  if />  n.  If  /  <  n,  then  let  v(s,  a)  be  an  arbitrary,  monic  polynomial  in  s 
with  coefficients  that  are  real  functions  of  a  such  that  at  each  aeT,  Sv(s,  a)  -  n  -f  and 
v  (s,  a)  has  only  negative  real  part  roots.  Then  we  can  write 

«)  nco  r 

df{s,  a)  "  d{s,  a)  ’ 
where  6df(s,  a)  -f>n. 

Step  3:  Assume  that  the  observability  index  of  a  minimal  parameterized  linear  realization  of 
Ta(s),  that  is,  the  row  index  of  Ta(s)  is  a  constant  v  for  aeT.  Let  de(s,  a)  be  a  monic  poly¬ 
nomial  in  s  with  coefficients  that  are  real  functions  of  a  such  that  at  each  aeT, 
8dc(s,  a)  A  r  >  v  - 1, 7 Tf(s,  a)/dc(s,  a)  is  proper  and  de(s,  a)  has  only  negative  real  part  roots. 
Set  w(s,  a)  ■  1  if  />  n  +  r,  and  if  /  <  n  +  r,  let  w(s,  a)  be  a  monic  polynomial  in  s  with 
coefficients  that  are  real  functions  of  a  such  that  at  each  a  €  I\  5h>  (r,  a)  *  n  +  r  -/<  r,  and 


w  (s,  a)  has  only  negative  real  part  roots. 

Step  4:  If  £[w>(j,  a)df(s,  a)]  <  n  +  r,  set  q(s,  a )  =  1.  Otherwise  use  polynomial  division  to 
write 

df(s,  a)  =  qi(s,  a)d(s,  a)  +  r2(s,  a)  ,  Srx(s,  a)  <  6d (s,  a) 

w(s,  a)qi(s,  a)  -  dc(s,  a)q(s,  a)  +  r2(s,  a) ,  Sr2(s,  a)  <  6dc(s,  a)  (11) 

In  this  division  process  it  can  happen  that  q\(s,  a)  will  have  coefficients  that  are  infinite  at 
particular  values  of  a  €  T.  A  similar  singularity  issue  can  arise  for  q  (s,  a),  and  also  the  possi¬ 
ble  occurrence  of  roots  at  s  =  0  in  q(s,  a)  complicates  the  subsequent  calculation  of  nonlinear 

compensators.  It  is  possible  that  these  difficulties  can  be  avoided  by  judicious  choice  of 

( 

dc(s,  a).  A  solution  in  any  case  is  to  increase  r  =  Sdc(s,  a)  to  avoid  the  divisions  in  (11), 
though  this  does  result  in  increased  compensator  dimension.  (Note  that  choosing  dc(s,  a), 
v  (s,  a)  and  w(s,  a)  to  be  independent  of  a  will  simplify  matters.) 

Step  5:  Let 

k(s,  a)  «  w(s,  a) nfa  a)  (12) 

and  determine  polynomials  l  (s,  a)  and  M(s,  a)  from 
w(s,  Q]df{s,  a)  -  dc(s ,  a)q(s,  a)d(s,  a) 

~  l(s,  a)d(s,  a)  +  M(s,  a )N(s,  a)  (13) 

Step  6:  Define  the  parameterized  linear  compensators  in  Figure  2  according  to 


& 

II 

F 

(14) 

c-(i> 

(15) 

M(s,  a) 
dc(s,  a) 

(16) 

It  is  straightforward  to  verify  that  the  parameterized  transfer  function  for  the  system  in 
Figure  2  is  7J(s).  It  should  be  noted  that  this  parameterized  linear  system  may  not  be 
minimal  in  the  sense  of  linear  systems.  However,  the  above  implementation  of  T^(s)  is  such 


that  the  system  will  be  internally  stable  so  long  as  (s)  is  stable. 

The  final  task  is  to  compute  a  nonlinear,  dynamic,  output  feedback  control  law  whose 
linearization  family  is  precisely  the  parameterized,  linear  control  law  shown  in  Figure  2.  Of 
course,  the  closed-loop,  parameterized  linear  system  in  Figure  2  described  by  7?  (s)  together 
with  {[w(a),  y(a)],  a€T}  is  a  linearization  family,  and  Ta(s)  with  {[u(a),  y(a)],  aeT}  is  a 
linearization  family.  Now  we  define  constant  operating  point  values  according  to 

■*«>  -  ■  «“>  - 

a  —(a) 

T(a)  -  /Cj(0)  %  da  ,  (17) 

*  I  3?W 

so  that  each  parameterized  linear  compensator  with  its  associated  constant  operating  point 
family  forms  a  linearization  family.  Using  (2)  and  (4)  it-  is  not  difficult  to  verify 
e(a)  =  r(a)  -  v(a)  explicitly,  that  is,  these  definitions  are  consistent  The  construction  of  non¬ 
linear  compensators  from  these  compensator  linearization  families  is  discussed  in  [5],  and  will 
be  illustrated  in  an  example. 

In  order  to  summarize  the  foregoing  development,  we  state  the  following  theorem. 

Theorem:  Suppose  we  are  given  a  model  (3)  and  a  nonlinear  system  with  linearization  family 
(1)  such  that,  for  aer,  where  T  is  a  sufficiently  small  neighborhood  of  o&,  d( 0,  a)  ^  0, 
N( 0,  a)  ^0,  dm{ 0,  a)  jfcO,  Nm(0,  a)  ^0,  and  (2),  (4)  hold.  Then  there  is  a  dynamic  output 
feedback  control  law  that  solves  the  linearized  model  matching  problem  if  for  a  e  r 

i)  Ta(s)  is  irreducible; 

ii)  the  row  index,  v,  of  Ta(s )  is  constant; 

iii)  N(s,  a)b  (s,  a)  =  a(s,  a)Nm(s,  a)  for  scalar  polynomials  a(s,  a)  and  b(s,  a), 
where  6a  (s,  a)  is  constant; 

iv)  6d  (s,  a)  -  6/N(s,  a)  <  6dm(s,  a)  -  5,Nm(s,  a),  i  =  1, . . . ,  p. 

Example:  The  dynamics  of  an  inverted  pendulum  on  a  cart  are  described  by 


*1 

*2 

x 2 

gsin(x!)  -  am£e|  sin^^cosCxi) -acos(x\)u 

4//3 -am/cos2(xi) 

y  -*» 


where  xj  is  the  angle  of  the  pendulum  from  the  vertical,  u  is  the  force  applied  to  the  cart, 
a  =  l/(m  +  Af),  m  is  the  pendulum  mass,  M  is  the  cart  mass,  2/  is  the  pendulum  length  and  g 
is  the  acceleration  due  to  gravity.  The  constant  operating  point  family  is  easily  parameterized 
by  the  first  component  of  the  state: 

u(ot)  *  ^-tan(a)  ,  x(a)  =  Q  ,  y(a)  =  a ,*  a€ T  ^  (-  j,  ~)  ,  oq  -  0 
The  transfer  function  of  the  linearized  system  is 


ar.(0  =  -I 


s2  -6(a) 


where 


6(a)- - L - ,  c(a) - — 

[4/  /3  -  am/cos2(a)]cos(a)  4//3  -  ami  cos2  (a) 

Note  that  ra(s)  is  irreducible  and  has  neither  pole  nor  zero  at  s  =  0  for  any  aeT  since 
41  /3  >  am/,  and  also  that  (2)  is  satisfied.  Let  the  model  be  specified  by 


"w'“ 

Obviously  (4)  is  satisfied,  and  the  precompensator  is  given  by 


7l(s)  = 


s2  -6(a) 


c(a)(s  +  l)2 


To  implement  T%,(s)  by  the  feedback  configuration  in  Figure  2,  we  apply  the  algorithm  as  fol¬ 
lows. 

Step  1:  Let  nfo,  a)  -  1  ,d/(s,  a)  =  c(a)(s  +  l)2. 

Step  2:  Since  Sd/(s,  a)  =  /  ■  2,  Sd(s,  a)  *  n  -  2,  let  v(j,  a)  ■  1,  and  set  df{s,  a)  -  df(s,  a), 
Jtf(s,  a)  =  n/s,  a),  with  /  *  /  -  2. 

Step  3:  Since  v  =  n  =  2,  let  r  *  1,  and  de(s,  a)  ■  s  +  1.  Then  /  <  n  +  r,  and  we  let 
w(s,  a)  =  s  +  1. 

Step  4:  Since  5[w(j,  a)df(s,  a)]  ■  n  +  r,  let  q  (s,  a)  *  1. 


Step  5:  Let  k(s,  a)  =  w(s,  a)rtf(s ,  a)  =  s  +  1,  and  solve  (12)  to  obtain 
/  (s,  a )  =  [c  (a)  -  l]s  +  3c  (a)  - 1  > 

M(s,  a)  =  [6(a)  +  3 >  +  36(a)  +  1 

Step  6:  The  parameterized  linear  compensators  in  Figure  2  are  easily  computed  to  be 


ci(s)  =  ci(5)  =  i 


From  (17),  r(a)  =  w(a),  e(a)  =  u(a),  and  v(a)  =  a  -  (g/a) tan(a).  Thus  the  nonlinear  com¬ 
pensators  associated  with  C„(s)  and  c£(s)  can  be  taken  as  unity  gains.  A  parameterized 

r(a)  -  -7a. 
v(a)  =  a  -  -^-tan(a) 

Finally,  it  is  easy  to  check  that  a  corresponding  nonlinear  compensator  is  described  by  [5] 

f  =  ,r.2y  +  %sm(y)-2aucos(y) 

41/3  -amlcos2(y) 

v  =  r-«  +  3v  +  gsinfrj-aHeosfr) 

4/  /3  -  ami  cos2  (y) 

Remark:  Linearized  model  matching  can  be  generalized  to  the  case  where  Ta(s )  and  T£  (j) 
have  poles  at  s  =  0  that  are  invariant  with  respect  to  a  6  T.  The  main  changes  to  be  made 
involve  a  more  complicated  characterization  of  the  associated  constant  operating  points  for 
the  linearization  families.  For  example,  suppose  Ta(s)  is  such  that  that 

d(s,  a)  =  (s”'1  +  di(a)sn'2  +  •  •  •  +  d„.i(a))s 

&d(s,a)s,  di(a)€R,  dn.1(a)  fO  for  a€T  (18) 

Then  we  must  have  u(a)  =  0  for  a  6  T.  To  characterize  the  associated  constant  operating 


linear  realization  for  the  compensator  C|(s)  is 
U  =  *  U  +  [ 2c  (a)  2(6  (a)  -  1)]  L 

vs  *  U  +  [c(a)  - 1  6(a)  +  3]|^J , 


point  output,  write  Ta(s)  as  the  concatenation  of 


Ta(s)  =  -  ^ — 

d(s,  a)(s  +  1) 

with  output  y  and  input  u,  and 


(19) 


Ta(s)  =  (20) 

with  input  u  and  output  u.  It  is  easy  to  see  that  Ta(s)  with  [u(a),  y(a)]  is  a  linearization  family 

-  A 

if  and  only  if  u(a)  can  be  found  such  that  Ta(s)  with  [u(a),  y(ar)],  and  Ta(s)  with  [u(a),  u(a)] 

A  #  5  . 

both  are  linearization  families.  Since  Ta(s)  with  u(a)  =  0  and  any  u(a)  is  a  linearization  fam¬ 
ily,  Ta(s)  with  {[u(a)  =  0,  y(a)],  a  e  T}  is  a  linearization  family  if  u(a)  can  be  found  such  that 


*L 

da 


(«)  -  r„(  0)|i(«), 


for  ttCf 


(21) 


Conclusions 

The  linearized  model  matching  problem  is  a  natural  extension  of  eigenvalue  placement 
problems  formulated  in  terms  of  extended  linearization  [1]  or  pseudo-compensation  [4]  in  that 
numerator  gains  and  zero  locations  of  the  linearized  closed-loop  system  are  explicitly 
addressed.  This  additional  capability  can  be  useful  when  acceptable  performance  cannot  be 
achieved  by  assignment  of  linearized-system  eigenvalues  alone. 

If  a  parameterized  compensator  can  be  found  such  that  the  model  can  be  taken  to  be 
independent  of  ct,  as  in  the  Example,  then  the  linearized  model  matching  problem  leads  to 
what  might  be  called  input-output  pseudolinearization.  However  this  bears  only  loose  resem¬ 
blance  to  the  original  notion  of  pseudo-linearization  in  that  variable  changes  are  not  involved, 
full  state  feedback  is  not  used,  and  the  dimension  and  eigenvalue  locations  of  the  closed-loop 
linearization  family  are  not  specified  a  priori 

One  disadvantage  of  the  frequency-domain  approach  used  here  is  that  extension  of  the 
results  to  the  multi-input  case,  though  clear  in  broad  outline,  is  difficult  in  detail.  For  exam¬ 
ple,  suppose  the  parameterized  linear  compensators  are  specified  in  terms  of  parameterized 
transfer  function  matrices.  In  the  course  of  constructing  the  nonlinear  compensators,  compu¬ 
tation  of  parameterized  linear  realizations  of  the  transfer  function  descriptions  can  lead  to 
singularities  with  respect  to  the  parameter  vector  a  that  are  not  easily  shown  to  be  removable. 


Also,  the  so-called  integrability  conditions  for  the  existence  of  a  nonlinear  control  law 

corresponding  to  a  parameterized  linear  control  law  are  nontrivial  in  the  multi-input  case. 
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Figure  1.  Linearized  plant  with  precompensator. 


